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Abstract 

We study properties of the zero and near-zero eigenmodes of the 
overlap Dirac operator in compact U{1) gauge theory. In the confine- 
ment phase the exact zero-modes are locaHzed as found by studying 
the values of the inverse participation ratio and other features. Non- 
zero-eigenmodes are less localized in the confinement phase. In the 
Coulomb phase no zero-modes are observed and the eigenmodes show 
no localization at all. 
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1 Motivation 



The Atiyah-Singer index theorem jTj relates the topological charge of the 
background gauge field configuration to the number of fermionic zero-modes 
of the Dirac operator. The theorem states a particular connection between 
the background gauge field and the fermionic fields. It is a theorem derived 
on differentiable manifolds, however. Quantum field configurations and in 
particular lattice fields are non-differentiable. Hence it is of interest to study 
on the lattice the topology of the gauge background as well as the properties 
of the zero-modes of the Dirac operator. 

A major progress in understanding the manifestation of the index theorem 
on the lattice was the realization that the eigenvectors of a 75-hermitian lat- 
tice Dirac operator with real eigenvalues should be interpreted as the lattice 
counterparts of the continuum zero-modes. This can be understood from the 
fact that only eigenvectors ip with real eigenvalues can have non-vanishing 
pseudoscalar matrix elements (ipl'^^l'ip), like the zero-modes in the continuum. 

Continuum QED in 4D in its usual non-compact realization for the gauge 
field has no topological charge. For the lattice formulation with compact 
representation of the gauge fields on the other hand the situation changes. 
Lattice QED exhibits a twofold phase structure: a physical one, containing 
the massless photon and a confining phase with properties similar to QCD. 
Some phenomena appearing in the confinement phase are: zero modes, mag- 
netic monopole-antimonopole pairs and the occurrence of a non-zero chiral 
condensate. 

Chiral symmetry breaking is a key feature of the theory of strong interac- 
tions, QCD. The order parameter is the chiral condensate (ipip)] it is related 
to the density of eigenvalues p{X) of the Dirac operator near the origin via 
the Banks-Casher relation 2J. As compact lattice QED shows a confining 
phase for certain values of the coupling constant one also expects jS] and 
indeed observes the appearance of a chiral condensate. 

For a study of zero-modes and chiral symmetry breaking it is advisable to 
work with chirally symmetric Dirac operators. On the lattice chiral symme- 
try is realizable in a locally violated form, expressed through the Ginsparg- 
Wilson relation |4j which in its original form reads: 

+ = 2aD75i?P . (1) 

where T> denotes the massless Dirac operator, a is the lattice spacing and 
R a local function of the gauge fields. In the continuum limit V anticom- 
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mutes with 75 thus showing chiral symmetry. Dirac operators satisfying the 
Ginsparg- Wilson relation preserve a lattice version of chiral symmetry j3] 
without fermion doubling. For a given local operator R there may be many 
Dirac operators. The eigenvalue spectrum of a Dirac operator for -R = | lies 
on a unit circle around the center 1/a and complex eigenvalues come in com- 
plex conjugate pairs. Only for real A the expectation value of 75 between the 
eigenstates (V'aItsIV'a) is non-zero. Therefore exact zero-modes have definite 
chirality. 

Although there are several realizations of the Dirac operator satisfying the 
Ginsparg- Wilson relation approximately, only the numerically costly overlap 
operator [U] is an exact realization. Within QCD there have been several 
studies on the density of near-zero-modes and the properties of exact zero- 
modes and their relationship to topological excitations (see e.g. |71IHII^ITU]). 

In compact lattice QED there have been studies with the overlap operator 
[TT] demonstrating the existence of zero-modes in the confined phase and 
relating the density of the near-zero-modes to the universal distributions from 
Random Matrix Theory. Here we want to continue the study of such modes 
with an emphasis on the chirality and the localization properties of exact 
zero-modes and possible correlations to topological structures like monopoles. 

2 Formalism 

The Wilson gauge action reads 

S[U]=f3 (1-cos^,,^,) , (2) 

X,fJ,>U 

where the gauge fields are represented by group elements U^^^ = exp^iO^^^) G 
f/(l) with 9^^^ G (— 7r,7r]. The plaquette angles are given by 6'^,.^^,^ = 9^.^^ + 
9x+fi,v — 9x+u,ii — 9x,u € (— 47r,47r). This compact realization of the gauge 
fields and the action introduces higher order self-interactions. The pure gauge 
theory has a confinement phase for small (3 and a Coulomb phase above a 
phase transition which for this action is close to /? ~ 1 and weakly 1st order 
O Cn] (for a recent study of the phase transition suggesting a new order 
parameter cf. Ref. One observes an abundance of monopoles in the 

confined phase. Monopoles in compact U{1) are defined following Ref. J3] 
from deficit angles of the plaquettes bordering 3-cubes, corresponding to links 
on the dual lattices. Due to current conservation the links form closed loops 
(better: networks ^Hl) on the dual lattice. 
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The massless overlap operator D^v may be written 



ov 



(3) 



where e denotes the operator sit^n-function and H = 75(3 — Hq) is some 
Hermitian Dirac operator, constructed from an arbitrary Dirac operator. It is 
convenient to use for Hq the usual Wilson Dirac operator with negative mass 
term and s is a real parameter which can be adjusted such as to optimize the 
convergence in the construction of e{H). If Hq is already an overlap operator 
then X>ov = Hq. 

For a U{1) gauge theory on manifolds with torus topology there is no 
topological charge in the sense of a non-vanishing Pontryagin index. So 
the Atiyah-Singer index theorem may be realized only trivially, i.e. allowing 
for a cancellation between the numbers of left-handed and right-handed zero- 
modes. However, on the lattice we have non-differentiable fields and therefore 
we may expect other violations as well. 

It has been demonstrated in QCD that for instanton configurations the 
Dirac operator shows a zero- mode localized in space-time [3 [TTl UHl E] • To 
our knowledge nothing is known about the localization properties of zero- 
modes in QED. 

In order to study the localization properties a simple set of gauge invariant 
quantities has been introduced. For an eigenvector iIj{x) of the lattice Dirac 
operator one defines a local density 



d 

where is an element of the Clifford algebra and the local sum runs over 
the Dirac indices. The eigenvectors are normalized to unit norm, i.e. 



For the cases of particular interest we abbreviate the scalar density p{x) = 
Po{x) (for Fo the unit matrix) and the chiral density p^i^x) (for 75). We also 
determined the vector densities p^ (for F = 7^), the axial vector densities 
P5^ (for F = 757^) and the tensor densities p^^ (for F = 7^71.^;,). 

Exact zero-modes are also eigenmodes of 75 with eigenvalues ±1. There- 
fore for such eigenmodes we expect even locally 




(4) 




(5) 



X 




(6) 
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The integral over the chiral density provides a measure for the amount of 
chiral symmetry breaking because it takes on its largest value for the zero- 
modes, namely ibl for normalized eigenvectors, and vanishes otherwise. 

The so-called inverse participation ratio (IPR) is introduced for further 
quantification of the localization (see for example Ref. [7j). We define it for 
both, p and p^. 

I = VY.p{x)\ h = VY.P^{^)\ (7) 

X X 

where V is the lattice volume. 

Assume that the density is evenly distributed on a subvolume Vf with 
p{x) = 1/Vf and vanishing elsewhere. Then we find / = V/Vf, with the 
limiting cases 1 = 1 for distribution over the whole lattice and I = V for 
localization at one point. It is essentially the inverse fraction of the volume 
contributing to the mode. This makes it an appropriate measure for the 
localization of eigenmodes. 

Let us try to understand why the inverse participation ratio may be 
considerably smaller than / for near-zero-modes than for exact zero-modes. 
The local Dirac sum may be written as a sum of two positive terms p{x) = 
where p5(x) = p+{x)—p-{x). Therefore > |p5(a:)|. Due to 
this property J5 is bounded for all eigenmodes: I > I5. For exact eigenmodes 
we have I = I5 due to (jHI). 

Consider now eigenmodes corresponding to eigenvalues far from the ori- 
gin: P5{x) will be much smaller than p{x) and fluctuate around zero. In this 
case I5 is expected to be significantly smaller than I. When the eigenvalues 
approach zero I5 increases and the ratio I5/I is expected to approach 1. 

We also determine different densities p„ for exact zero-modes of chirality 
+ 1 or -1. This amounts to checking relations like e.g. 

(^ItmTsV') = M^M , (8) 

which we could verify. 

3 Simulation and results 

The background gauge fields were obtained by using the Metropohs and over- 
relaxation updating algorithm. All configurations are well decorrelated, sep- 
arated by 5000 updating sweeps. We analyzed 400 configurations on the 4^, 
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Figure 1: The plots show the normahzed histogram (probabihty distribution) 
for the topological charge v for the configurations at /5 = 0.9 (left-hand plot) 
and /3 = 0.99 (right-hand plot). 



6^ and 8^ lattices and 100 on 12^ lattices. For the computation of the eigen- 
values and eigenvectors the so-called implicitly restarted Arnoldi method 
[inj was used. The overlap operator was computed by an appropriate se- 
ries expansion. Only the smallest 10 eigenvalues and their eigenvectors were 
calculated. Real eigenmodes at zero have exact values ('?/'|75'?/') = ±1. 

In the confined phase many exact zero-modes were observed whereas no 
zero-modes appeared in the Coulomb phase. Zero modes are chiral eigen- 
modes and are related to the topological charge via the Atiyah-Singer index 
theorem. In our discussion we therefore identify the number of zero-modes 
with positive or negative chirality in a configuration with the positive or 
negative topological charge v. We never observed configurations which have 
zero-modes of different chirality. To our knowledge there is no formal theo- 
rem explaining this feature in 4D (there exists such a vanishing theorem in 

2D [2ni). 

In Fig. ^ we show the number v of exact zero-modes with chirality ±1 
for two different values of /?. We observe obvious volume dependence. The 
number of zero modes as well as the average percentage of zero-mode configu- 
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Table 1: Summary of configuration number with topological charge for (3 = 
0.99 and the corresponding topological susceptibihty. 

rations compared to all configurations grows with the volume. The integrated 
number of zero-modes (i.e. adding up all |z/|) rizeros/nconfigs reaches values 
(9(0.9) for the large lattices. This can be observed for both /3 = 0.99 and 
P = 0.9. On 4^ lattices at /3 = 0.99 we found only 3 such modes on all 400 
configurations (compared to e.g. 24 zero- modes on 6^). This small num- 
ber can possibly be attributed to the situation that the pseudocritical point 
(i.e. the peak position of the specific heat and other cumulants) on smaller 
lattices effectively occurs already at smaller /3-values and thus we may have 
been already in or at least closer to the Coulomb phase for that lattice size. 
For 6^ the position of the pseudocritical point is near 1.002 whereas for 12^ 
it is near 1.010, cf. Ref. \n\ . 

Table n summarizes the results for j3 = 0.99 together with the topological 
susceptibility. Its behavior is compatible with approaching a constant for 
large volumes. 

Following the definitions from Refs. ^31^ we determined the monopoles 
(which are closed loops on the dual lattices). We find no correlation between 
the monopole density and the number of exact zero-modes on the corre- 
sponding configurations. For (3 = 0.9 the monopole density is 0.19(1) for 
both lattices sizes 6^ and 8^, for (3 = 0.99 it varies from 0.118 to 0.121 for |z/| 
ranging from to 3 on the largest lattice size 12^. For (3 > (3c the density 
decreases exponentially. We see no correlation between the number of zero- 
modes and the integrated monopole loop lengths at /3 = 0.9 and 0.99; this 
observation confirms jTT]. 

The scaling of the density of small near-zero-modes in the sector of dif- 
ferent zero mode number has been compared with Random Matrix Theory 
expectations in Ref. jTT]. The universality class was identified to be the 
unitary ensemble. We also studied the densities for our data and agree with 
these findings. Due to our choice of a different offset parameter value s in 
the Wilson action entering the overlap operator construction, the condensate 
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Figure 2: We compare the densities of the smallest eigenvector for (3 = 0.9 
(in the i/ = sector) and (3 = 1.03 and 1.1 for lattices size 8^. 

value (and with it the scaling parameter) has to be renormalized multiplica- 
tively [22] for the comparison. 

Fig. |21compares the densities for the confinement region with those for the 
Coulomb region, exhibiting a different behavior. In the confinement region 
one expects for z/ = a distribution following the chiral unitary ensemble of 
Random Matrix Theory [2H] and we overlay a fit to that functional behavior 

= £ exp (-zV4) . (9) 

In the Coulomb region the results for j3 = 1.03 and (3 = 1.1 are compatible 
with each other indicating that there is no gap, different to the observed 
behavior in the deconfined phase in QCD. Our statistics is not sufficient to 
exclude a small gap, though. 

For the zero- modes we observe p{x) = ±p^{x) as expected for exact GW- 
fermions and thus also the IPRs agree. For chiral eigenmodes we expect local 
agreement of those densities which correspond to relations for the gamma 
matrices like 7172 = —737475, i-e. ^12(2;) = — P34(x). Indeed the tensor 
densities for the chiral modes do obey this relation within their group. The 
local densities for vector and axial vector both vanish. For non-zero-modes 
these relations are not valid. 

In general the integrated tensor densities are comparatively small in ab- 
solute magnitude (0.1... 0.01) for zero-modes. For near- zero-modes all in- 
tegrated vector, and tensor densities are small (0.1 . . .0.01). The integrated 
pseudoscalar density was exactly ±1 as expected for exact zero- modes. 
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Figure 3: The graphs show the mean inverse participation ratio for the near- zero- 
modes (in the = sector) at P = 0.99 and /3 = 1.03 (lattice size 8^): mean 
values of IPR (l.h.s.) and IPR- histograms (r.h.s.). The exact zero modes in the 
= 1 sector have a wider distribution (shaded area) extending to higher values; 
the (/) and s.d. for those is indicated by the fihed circle (arrow) and error bars in 
the upper l.h.s. plot. 

The inverse participation ratio should allow some information on the 
space-time localization of the eigenmodes. In Fig. Olwe summarize this quan- 
tity for the modes for non-zero eigenvalues, binned according to Im(A) and 
also show the IPR histograms. In the confinement region we find no increase 
of the localization for the lowest bins. This observation differs from results in 
QCD studies f7j, where also the near- zero modes show stronger localization. 

In the Coulomb region no modes at very small Im(A) were observed and 
the low- lying modes show small IPR, close to the minimum of 1. 

Let us now concentrate on the exact zero-modes in the confining phase. 
These are well localized objects with comparatively large individual IPR 
values up to 4.2. For fixed (3 the (J) show little if any volume dependence 
with values of 1.92(29) for L = 8 and 2.03(69) for L = 12 (this independence 
on the volume was also observed in recent QCD studies These values 

are, however, clearly larger than the mean values for the small non-zero 
modes. The difference is mainly due to an IPR distribution extending to 
much higher individual IPR values, as seen in Fig. El The values of IPR for 
the lowest non-zero modes for the sector with no or that with one exact zero 
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Figure 4: The density shape for zero-mode eigenvector of a Monte Carlo 
generated gauge configuration (at (3 = 0.99, lattice size 8^). The 3D cuts 
are for time slices = 1 . . . 4 (first row) and = 5 . . . 8 (second row) . This 
zero-modes has / = 2.75. 

mode agree within small errors. 

In an attempt to visualize the geometrical shape of these localized zero- 
modes we studied 3D cuts of the density distribution. Fig. Ogives an example 
of the observed structure. We plot the surfaces of constant density p{x) for 
an exact zero-mode for the eight 3D time slices (nj = 1 . . . 8) of the 8^ lattice. 
The shape describes a tubular structure in 4D. However, most of the observed 
shapes are less clear and often made up from disconnected pieces. Contrary 
to QCD PB] no clear instanton-like picture (i.e. 4D "blobs") evolves. The 
structure we observe for the zero-mode eigenvectors are sometimes ball-like 
in 4D, sometimes tube-like, even tubes closed in some direction due to the 
periodic boundary conditions. Although a relation to the concept of closed 
monopole loops (cf. Ref. ^T) may be conjectured, no stringent conclusion 
can be drawn. 

This evidence supports the notion, that the zero-modes are not instanton- 
like lumps, but extended and periodically closed in one or more directions. In 
a recent study of the zero-modes in QCD it was argued, that the scaling 
of (/) with the lattice spacing may be used to identify the co-dimension of 
these objects. Due to lack of results at several different lattice spacing we 
cannot attempt this here. In Ref. values of (/) ~ 2.8 were observed, i.e. 
larger than our values (!) ~ 2, indicating even stronger localization. 

In order to gain further insight, we determined for some gauge config- 
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Figure 5: The density shape for zero- mode eigenvector of the constructed 
gauge configuration (lattice size 8^). The 3D cuts are for shces rij = 1 . . .4 
(first row) and rit = 5 . . .8 (second row). The inverse participation ratio for 
that mode is / = 1.85. 

urations the eigenmodes of the Dirac operator applying also anti-periodic 
boundary conditions in the time direction by multiplying all temporal gauge 
links of a timeslice with -1. This leaves the gauge action invariant (although it 
is not a gauge transformation on the finite lattice). The number of monopoles 
and Dirac plaquettes does not change, either. 

However, the characteristic polynomial coefficients of the eigenvalue equa- 
tion for the Dirac matrix involve traces over closed loops and thus may differ 
due to the periodically closed loops. The Dirac matrix was diagonalized again 
for the such modified gauge configuration and the eigenvalues were compared 
with the ones without an added phase: Quite often the zero-modes disap- 
peared. Where the zero-mode survived the transformation, the 3-D density of 
the new eigenvectors had a completely different structure. This confirms the 
suspicion that for U{1) gauge theory the zero modes are related to lumps ex- 
tending in at least one direction. Changing the boundary conditions in QCD 
sometimes also repositions topological objects (see e.g. |25J and references 
therein) and in a recent analytical study it was reported [2S] that zero modes 
with constant curvature may be constructed in the U{1) subgroup. These 
modes may be "switched" on and off by changing the boundary conditions. 
This is in agreement with our observation. 

In 2D the Pontryagin index is constructed from / da^^F^y. On the lat- 
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tice for compact U{1) gauge theory and with torus topology various exphcit 
configurations with such topological charge have been explicitly constructed 
as well as observed in Monte Carlo simulations (e.g. in Ref. |2Z1)- In 4D 
the topological charge is proportional to / d'^xFfj^iyF^^ and a possible method 
123 to construct objects with non- vanishing topological charge is e.g. to 
combine two-dimensional submanifolds, both with topological charge in the 
2D definition. Intersections of 2D surfaces in 4D may be points as well as 
more complicated, e.g. line-like objects. In Ref. JT] such a configuration - 
as suggested in Ref. [21] has been studied and identifies as leading indeed to 
exact zero-modes. For unit topological charge (lattice size L"^, uj = 27r/L^) it 
has the form 

Ui{x) = exp (207x2), (10) 
U2ix) = 1 for X2 = 1 , (11) 

exp(— ztuLxi) for a;2 = , (12) 

and for f/3 and f/4 equivalent, with 0:3 and replacing xi and X2- Fig. El 
shows the geometric shape of the eigenvector density of this mode and we 
observe tubular structures. Note, however, that the tubes are in the 3D 
intersections, i.e. extend like planes in 4D. 

4 Conclusion 

We have studied zero-mode properties of the overlap Dirac operator for f/(l) 
gauge theory in the confined phase. We find no correlation between the 
occurrence of these modes and the monopole density confirming [11]. 

We do observe individual zero-modes with definite chirality. This would 
contradict the Atiyah-Singer index theorem if applied naively for lattice sys- 
tems. However, for U{1) lattice gauge theory for compact gauge fields the 
situation changes. In particular the confinement phase of this system has 
no analytic continuation to the weak coupling continuum limit. In the con- 
finement phase one can explictly construct configurations with topological 
charge and zero-modes [21], essentially by combining two 2D sub-manifolds 
with non-vanishing 1st Chern number. 

In the Coulomb phase there are no zero-modes and also the near-zero- 
modes are suppressed. There we find no pronounced localization signal. In 
the confinement phase the exact zero-modes are definitely more localized than 
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other eigenmodes (as exhibited by the inverse participation ratio). The IPR 
of the nearby non-zero modes does not seem to depend on the Im(A). The 
geometrical structure of the zero eigenmodes has no clear signature singling 
out 4D blobs - in contrast to e.g. the case of instantons in QCD, which 
have been observed in distorted hyperspherical shapes. All evidence points 
to tubular or even planar structures supporting the density function p{x) for 
a zero-mode. 
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